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Abstract 

Within Tsallis' nonextensive statistics, a model is elaborated to address self-similar 
time series as a thermodynamic system. Thermodynamic-type characteristics rel- 
evant to temperature, pressure, entropy, internal and free energies are introduced 
and tested. Stability conditions of time series analysis are discussed in details on 
the basis of Van der Waals model. 
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1 Introduction 



Time series analysis allows to elaborate and verify macroscopic models of complex systems 
behaviour on the basis of data analysis. The latter is known to be focused on numerical 
calculations of correlation sum for delay vectors that allow to find relevant magnitudes of 
fractal dimension and entropy [1] . Being traditionally a branch of the theory of statistics, time 
series analysis is based on the class of models of harmonic oscillator, which are relevant to the 
simplest case of the Gaussian random process [2] . But it is well known that real time series is 
rather the Levy stable processes related to self-similar system, than the Gaussian ones being 
very special case [3]. Therefore, it is important to study self-similar time series. In this Letter 
we present a method based on Tsallis' statistics [4] , which is of relevance for self-similar system 
analysis [5] - [7]. 

The paper is organized as follows. Section 2 is devoted to elaboration of a model, which 
permits to address a self-similar time series as a nonextensive thermodynamic system. Section 
3 is based on calculations of both entropy and internal energy of the time series. As a result, 
thermodynamic-type characteristics of the time series such as temperature and entropy, volume 
and pressure, internal and free energies are introduced. Their testing for the model of ideal gas 
is shown to be basis for statistics of self-similar time series. Section 4 is devoted to discussion of 
the physical meaning of the results obtained. Stability conditions of time series mimicked as a 
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on. Several equalities needed in quoting are placed in Appendix. 

2 Presenting time series as nonextensive thermodynamic system 



Let us consider D-dimensional time series x(i n ) related to the set {x n } of consequent values 
x n = x(t n ) of principle variable x(i) taken at discrete time instants t n = nr that we obtain as 
result of dividing a whole time series length T = Nr by N equal intervals r. It is obviously 
to be relevant to the time series x(t n ), the set {x n } should be supplemented by conjugated 
set {v„} of velocities, which show rates of x n - variation with the time jumping. In the simplest 
case of Markovian consequence, one has v n = (x n — x n _i)/r. For more complicated series with 
m-step memory, the velocity magnitude is defined as follows: 



v n = 



lib ■ -t 



$i(m, n) = X(w " m)+i x ("-™)+(*-D , (i) 



The paradigm of our approach is to address the time series as a physical system defined by an 
effective Hamiltonian H = 7i{x„,v„}, on whose basis statistical characteristics of this series 
could be found. If one proposes that series terms x n related to different n are not connected, 
the effective Hamiltonian is additive: 

N 

H=^2 £ n, £n = (x n ,V n ). (2) 
n=l 

Physically, this means that the series under consideration is relevant to an ideal gas comprising 
of iV identical particles with energies e n . Further, we suppose different terms of time series to 
be statistically identical, so that effective particle energy does not depend on coordinate x n : 
s(x n , v n ) =3- s(v n ). Moreover, since this energy does not vary with inversion of the coordinate 
jumps x n — x n _i, the function e(v„) should be even. We use the simplest square form 

?n = ^1, (3) 

which is reduced to the usual kinetic energy for a particle with mass 1. With switching on an 
external force F = const, particle energy (3) becomes as follows: 

£n = - Fx n . (4) 

Finally, when time series has a microscopic memory, dimension of the delay vectors 5i(m,n) 
in the definition (1) needs taking m > 1. Moreover, if time series terms x m , x n with m ^ n are 
clustered, the Hamiltonian becomes relevant to a non-ideal gas with interaction w mn , m ^ n: 

w = ^IX + ^E«w (5) 

n=l m^n 
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(0) q = J2 0{x„,v n }P g {x n ,v n } 

{x n ,v n } 



(6) 



defined by relevant functional 0{x n ,v n }. Hereafter, the principle role is devoted to the sum 
over states, which reads: 



E 

{x n ,v„} 



1 

ivT 



dx n dv n 



n 

n=l 
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N(tA) d 



n / d y« / du « =► M' 1 II / d y™ / du - ( 7 ) 

-, |dy n = |du n = l.(8) 



yn = Y' u ^ 



X 



Here, A is effective Planck constant that determines a minimal volume of phase space per a 
particle related to a term of the time series, the factorial takes into account statistical identity 
of these terms, the factor (X 2 /t) dn is caused by change of initial variables x n , v n into y n , u n 
rescaled with respect to macroscopic length X. 

At a choice of the probability distribution P g {x n , v n } =>- V q {y n , u n }, we are stated on a self- 
similarity of the time series meaning a power-law form of the functional V q {y n , u n } instead of 
usual exponential one type of Gaussian. We adopt known Tsallis' escort distribution [4] 



V q {y n , u n } = < 



at (l-g) ^gP >l, 



1 _ M - a) Myn,Un}-g 
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otherwise. 



(9) 



Here partition function is defined by condition 



Z = 



N 

= n / 

n=l J 



l-(l-9) 



7i{y n ,u n } - £ 



9 
1-9 



dy n du n , 



(10) 



where < q < 1 is a parameter of nonextensivity, T is energy scale, internal energy E and 
normalization parameter (1) are determined by equalities 



N r 

E = M 1 Y[ / H{y n , u n }V q {y n , u n }dy n du n , 

n=l J 



N f 

II / CPqiyn, u n })" dy n du r , 

n=l n 



(11) 



To check the statistical scheme proposed let us address firstly trivial case of time series x„ = 
const. Here, the particle energy e is a constant as well, so that the Hamiltonian is 7i = Ne. 
The partition function Z = A/" -1 and the normalization parameter (1) = J\f~^~ q ^ are given 
by inverted normalization factor (8), whereas the internal energy E = Ne is reduced to the 
Hamiltonian. Then, the entropy H = — In M obtained according to definition (59) given in 
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we find effective Planck constant: 




NJ r 



(12) 



If D-dimensional domain of the coordinate x n - variation supposes to be governed by Levy-type 
law 



with microscopic constant x and dynamic exponent z, then one obtains the following scaling 
relation for the phase space volume per a term of the time series: 



In the case of Gaussian scattering, when z = 2, this volume does not depend on number iV of 
the time series terms, i. e., one obtains a constant. 



3 Non-extensive thermodynamics of time series as an ideal gas 

Stating on the basis of known results for thermodynamic quantities of nonextensive ideal 
gas (see Appendix) we are in position to examine pseudo-thermodynamic properties of time 
series. The quoted equalities (49) - (52) derive to the starting expression for the normalization 
quantity 



x D = x D m 



(13) 




D 



(14) 




l-o 



> 



(15) 



Here the notations are introduced (cf. Eqs. (51)) 




a = 




(16) 



where we use T-function. In the limits iV>l, (l-g)|<l, when 




(17) 



one obtains for the entropy (59): 




(18) 
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^ = iV^ln(§), G = (2vreT )^g)^, (19) 
if the dynamic exponent is determined as 



(20) 



l-a 

Respectively, the internal energy (50) and the normalization parameter (15) read: 
DN /G\*d-i) ,., /G\U 1 -i) 



E = — t »- ®<={n (21) 



Following [8], let us introduce the temperature 



/G\i( 1 - 1 *) z-1 

t = W* t ° = {n) To > a = — W 

where the second equality takes into account the last relation (21). This definition guarantees 
the equipartition law 



E = CT, C = cN, c= ^, (23) 
where the quantity 



C = ^ (24) 

is the specific heat. It is easily to convince that equations (19) - (22) arrive at standard 
thermodynamic relation 



— = - (25) 

Above used treatment is relevant to the method [9] addressed to a fixed value of the internal 
energy E. In alternative case when the principle state parameter is the temperature T, we 
should pass to the conjugate formalism [10]. Here, standard definition 

F = E-TH (26) 
of the free energy arrives at the dependence 
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F = - CTln UJ- (27) 

Then, the thermodynamic identity 
dF 

W - -B (28) 
yields the relation 

tf = Cln(£) (29) 

that plays a role of a heat equation of state. It arrives at the usual definition of the specific 
heat (cf. Eq. (24)) 

a tt 

C = T-. (30) 
Let us introduce now specific entropy per unit time 

h=(Dr)- 1 ^=T- 1 H 1 -r (31) 
to be determined by a minimal entropy 

#i = z In y/2neT - z In (^j (32) 
and a redundancy 

r = ( J Dr)- 1 ln(eiV). (33) 
Dependencies on the scale x 

z ( x \ 

h(x) = const In ( — J , r(x) = const (34) 

notice that the system behaves in a stochastic manner [11]. 
Let an effective pressure be defined as 

p = -r- (35) 

to measure specific entropy variation with respect to the time series scale. Then, we arrive at 
a mechanic equation of state 



being additional to the relation (29) of entropy to temperature. According to Eq. (36), defini- 
tion of the compressibility 



(37) 




dp J 

shows that it is reduced to the dynamic exponent (20). Respectively, stability conditions 

C > 0, k > (38) 
reduces to natural restriction 

z>\. (39) 

It is principally important that the "pressure" (35) is introduced as derivative of the specific 
entropy h with respect to the microscopic scale x. This means a microscopic nature of so 
defined "pressure", which related to a susceptibility (more exactly, the compressibility (37)) 
of the time series with respect to a choice of the microscopic scale x. 

4 Discussion 



We have addressed above the simplest model of the ideal gas, which have allowed us to examine 
analytically a self-similar time series in standard statistical manner. Characteristic peculiarity 
of related equalities (12) - (14), (18) - (21), (27), (29), (32) - (34), (36), (38), (39) is a scale 
invariance with respect to variation of the nonextensivity parameter 1 — q, which is contained 
everywhere as combination (see Eq. (16)) 

DN 

a = (l-q)C C=— (40) 

with free magnitude of parameter a related to the dynamic exponent z (see Eq. (20)). This 
invariance is clear to be caused by self-similarity of the system under consideration. For a given 
time series, the value of exponent z is fixed if it is addressed to a monofractal manifold and 
takes a closed set of magnitudes z in the case of self-similar system relevant to a multifractal. 
In real time series this invariance can be broken, so that a dependence on the parameter of 
nonextensivity could appear. However, above introduced set of pseudo-thermodynamic char- 
acteristics of time series is kept as applicable and accustomed thermodynamic relations (23), 
(24), (28), (30), (35), (37) can be applied to analysis of arbitrary time series. 

Let us focus now on effect of external field and particle interaction, which switching on is 
expressed by equalities (4), (5). Related entropy additions (54), (56) arrive at total value of 
the specific heat (30) in the following form: 

C,^0- C -{FX-^ w ), (4!) 
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T c = FxN^ - -J^ wx ~° Nl ~~^ ( 42 ) 

the system becomes nonstable due to external force F and interparticle attraction w < 0. If 
the number _V is much less than critical value _V C defined as 

^ 4)4-. __ £Z__M X -( D+ .), (43) 

main contribution gets interparticle attraction w < 0, in opposite case N ^> N c - external 
force F. With accounting temperature definition (22) and notation (19), this instability means 
that a microscopic time interval r should be more than a critical magnitude r c . In the case 
_V < _V C , we find 



r c = (^> l ) 2(Z " 1) x i-_ T ^ TjV A(iH). (44) 
In opposite case _V 3> _V C , one obtains 
i 

J?2(z-1) 2z-l 2z-l 

._ X 2!j=T) N 2D ^- 1 ) . (45) 



o 

On the other hand, the entropy additions (54), (56) accompanied by Eqs. (31), (35), (36) 
arrive at the compressibility (37) in the form 

«-' = *-(I> + l)^. (» + £). (46) 

Hence, the system becomes nonstable if microscopic scale x is less than a critical magnitude 
determined as 

x? = (D + (v + (47) 

In the limit of very large time intervals r, when the temperature is T ^> w/v, one finds the 
simple estimation 

x c ~ eN^ 1 '*). (48) 

Thus, at z = 1 critical value x c is reduced to a scale of the time series resolution e, being a 
particle radius within the gas model. 

Above studied instabilities break the conditions (38), which guarantee predictability of the 
time series analysis. Our main result is that the predictability becomes inpossible if microscopic 
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N of time series terms. 



To conclude, we point out that main advantage of the approach proposed is a possibility of its 
application to numerical analysis of real time series. The basis of this calculations is stated on 
expressions (1) - (6), (9) - (11), (22), (23), (24), (26), (28), (30), (31), (35), (37), (59). This 
work is in progress. 
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Appendix: Non-extensive thermodynamics of thermodynamic system 



Calculation of main thermodynamic quantities of nonextensive ideal gas arrives at the Euler 
T-function to be reduced to the following expressions for partition function and internal energy 
[8] -[10] 



V N l(q) 
N\ 



Hl) q 
1-q 



DN 

2 



! + (!-<?) 



DN 



E 



1 i DN 
2 



(49) 



E 



DN V N -i{q)T<„ 
~2 N~\ 



1-q 



DN 
2 



1 + (!-<?) 



DN 



1 -l 



E 



1 + 



1-9 1 2 



Z- q . (50) 



Here, D-dimensional gas in volume V = X D is addressed and the notations are introduced 



(51) 



v • r(^ + ^)' 

where m is particle mass, h is Planck constant and T is energy scale. Normalization quantity 



follows from the normalization condition of the Tsallis' distribution function [4]. 



(52) 



Switching on an external force F = const causes the second term in Hamiltonian (4) to arrive 
at the following factor in Eq. (49): 



Z ext = exp 



DN 



/FX 
\~2T~ 



smh (g 



DN 



(53) 
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/FX \ 



sinh 



(54) 



Within a supposition F > 0, a homogeneous external field arrives at monotonous decrease of 
the entropy H ext (T) with the temperature growth. At T > FX contribution of the second 
term is neglecting, in the limit T — > it tends to the first one. 

According to recent work [12], cluster expansion of the particle interaction in Hamiltonian (5) 
results in additional factor in partition function (49): 



N 2 / W\ „ 7 , , n , , „ 2ti d / 2 S 



Zint = 1 - 7^7 [V + — ; W = 



= S D J w(x)x D - 1 dx, S D ee ; v = ^e D , (55) 



where e is effective radius of the particle core. Relevant entropy addition 



aN 2 ( w\ 

H int = -—lv + -\ (56) 



21/ V T 

monotonously decreases with growth of the particle volume v. Growth of the temperature T 
causes the entropy decrease for attractive interaction w < and its increase in the case of 
repelling one (w > 0). 

Finally, we arrive at useful relations between physically defined entropy H and related Tsallis' 
magnitude H q . One has symbolic equation 

DN 

H = a\n[exp q (H q )], H q = lnjexp(#/a)]; a = (l-q) — , (57) 

that are alternations of the usual functions logarithm ln(x) and exponential exp(x) with cor- 
responding Tsallis' generalizations [4] 

X 1 ~ q — 1 _!_ 

\a. q (x) = - - - , expire) = [1 + (1 - q)xy-i. (58) 
In explicit form, relations (57) are appeared as 

H = alnZ, H q = ^ g ~ - . (59) 
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